Indirect magnetic interaction mediated by spin dimer in Cu2Fe2Ge40i3 
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Cu2Fe2Ge40i3 is a bicomponent compound that consists of Cu dimers and Fe chains with separate 
energy scale. By inelastic neutron scattering technique with high-energy resolution we observed the 
indirect Fe - Fe exchange coupling by way of the Cu dimers. The obtained parameters of the effective 
indirect interaction and related superexchange interactions are consistent with those estimated semi- 
statically. The consistency reveals that the Cu dimers play the role of nonmagnetic media in the 
indirect magnetic interaction. 

PACS numbers: 75.10.Jm, 75.25. +z, 75.50.Ee 
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Some types of magnetic interaction in solid are indi- 
rect ones by way of nonmagnetic media. The superex- 
change interaction among the most common type in in- 
sulating metal oxides is derived by perturbative treat- 
ment of electron transfer between nonmagnetic anions 
and magnetic ionsi. Another type is found in RKKY in- 
teraction in magnetic alloys^ where free electrons in the 
nonmagnetic metal transfer exchange integral between 
magnetic atoms. Anions or nonmagnetic metal behaves 
as the media for the indirect interaction. As the exotic 
types spin-singlet entities such as spin dimers, ladders, 
spin integer chains^ can play similar role. The theoret- 
ical treatment was considered in nonmagnetic impurity 
doped S = 1 chainsii^, S = 1/2 ladders^i, or 5 = 1/2 
dimerized chains^i^. Around the impurity the induced 
staggered spins generate an effective spin of S' = 1/2 in 
the sea of nonmagnetic state. The low energy excita- 
tion are described by the effective spins that indirectly 
interact with the similar formula to the RKKY interac- 
tion. This scenario is experimentally realized in doped 
spin gap compounds such as CuGeOa SrCu203^^, 
and PbNi2V20g^^. However, randomness of the impu- 
rity distribution makes the system rather complex. For 
the basic understanding a bicomponent spin compound 
with uniformly distributed magnetic ions is necessary. 

Recently Fe-doped CuGeOs with uniform Fe dis- 
tribution was identified in insulating metal oxide 
Cu2Fe2Ge40i3; Cu„_2Fe2Ge„03„+i with n = 4 in a 
hypothetical series of homologous structures with n — 
3,5,6, ...(X>i^. Here the compound with n = oo corre- 
sponds to well known spin-Peierls cuprite CuGeOs — 
where uniform ID chains of Cu^"*" ions run in the c direc- 
tion as is illustrated in Fig.[T](a). In Cu2Fe2Ge40i3 50% 
of Cu^"*" {S = 1/2) ions are periodically replaced by Fe'^+ 
{S = 5/2) and the uniform chains become distorted al- 
ternating chains of S* = 1/2 and S — 5/2 dimers in Fig.[T] 
(b). Consequently the lattice unit is quadrupled along 
the chain direction and the c-axis in CuGe03 changes to 
the a axis in Cu2Fe2Ge40i3. Bulk susceptibility and neu- 
tron inelastic scattering (NIS) measurements indicated 
that the Jcu-Fc are wea k^^'^^ and, instead, the strength 



of the interdimer interaction of Fe ions (J^) is the same as 
that of the intradimer one (J(,) within the experimental 
error. Hence Cu2Fe2Ge40i3 is a weakly coupled quan- 
tum dimers of Cu ions and classical uniform chains of Fe 
ions in the b direction. In this geometry Cu dimers seem 
to behave as nonmagnetic media for effective exchange 
interaction (Jcs) of neighboring Fe ions as is described 
in Fig. [1] (c). A set of exchange parameters were deter- 
mined by NIS and neutron diffraction (ND) in the previ- 
ous study as is summarized in the first line of Table [D In 
spite of the rigorous experiments the direct measurement 
of Jeff was not performed because of the twinning of the 
available crystal. The values in the square bracket is the 
reference estimates by a semi-static method. 

In this paper we performed cold neutron inelastic scat- 
tering experiment with high-energy resolution on the 
twinned but high-quality crystal. Thorough elaborate 
data analysis we successfully obtained the low energy ex- 
citation in the a* direction where the Jeff makes main 
contribution. We analyze whole dispersion in combina- 
tion with previously reported data^^ by both the con- 
ventional spin wave (SW) and the effective SW model 
based on the energy scale separation. We found the lat- 
ter model more appropriate in quantitative level. The 
obtained consistency between the Jeff's estimated by the 
dynamics measurement and the semi-static method re- 
veals that the Cu dimers play the role of the nonmagnetic 
media in the indirect interaction between Fe spins. 

The high quality crystal with the dimension of 7 x 5 x 40 
mm'^ was obtained by floating zone method. NIS exper- 
iment was performed on LTAS spectrometer installed in 
JRR-3 in JAEA. Open-80'-80'-open with soUar coUima- 
tion and fixed final energy of Ej — 3.5 meV were em- 
ployed. Closed He cycle refrigerator was used to achieve 
T = 3.0 K. The crystal structure of Cu2Fe2Ge40i3 is 
monoclinic. The obtained crystal was twinned so that 
a* and b* axes are shared as is illustrated in Fig. [1] (d) . 
Here the reciprocal spaces for crystal A and B are de- 
fined as a*b*c* (dashed grid) and a'*b'*c* (dashed dot- 
ted grid) coordinates, respectively. In Cu2Fe2Ge40i3 
a reciprocal vector {h k I) in crystal A corresponds to 
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FIG. 1: (Color online) (a) Schematic structure of 
CuGeOs projected onto c - b (left panel) and c - a (right 
panel) planes, (b) Schematic structure of Cu2Fe2Ge40i3. 
Lattice unit cell and magnetic unit cell are shown by solid 
and dashed lines, respectively. Main exchange path is crank- 
shaft Fe chains along b direction and Cu dimers. (c) Indirect 
Fe - Fe coupling by way of Cu dimers. (d) The structure of 
the twinned crystal of Cu2Fe2Ge40i3 projected onto a* - c* 
plane. Dot-dashed and dashed grids are for crystal A and B, 
respectively, (e) Schematic structure of Cu2Fe2Ge40i3 based 
on the effective SW model (see text). 



{-{h + 0.53/) - fc in crystal B. The a*b*c* coordinate 
will be used hereafter. 

A series of constant q scans were performed on the tra- 
jectory of the thick solid line in Fig. [1] (d) . Peak profiles 
of the typical scans are shown in Fig. [2] (a) . At antiferro- 
magnetic zone center of h = 2.0, well-defined spin wave 
excitation with an anisotropy gap of about 1 meV was 
observed in crystal A. The peak is wider than experimen- 
tal resolution and seems to have a structure probably due 
to xy two ion anisotropy. The shape of the scan profiles 
are well explained by considering resolution function in 
the Cooper-Nathans approximation^- that is shown by 
solid curves. AX fiio ^ 2.4 meV, broad peak from crystal 
B is detected. These two peaks move with the change 
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FIG. 2: (a) Constant q scans at (/i 1 — 0.5). The dashed- 
dotted and dashed curves indicate excitations from crystal A 
and B, respectively. Solid curves are fitting curves by resolu- 
tion convoluted function, (b) Dispersion curves of crystal A 
(dashed-dot) and crystal B (dashed). 



of h on the identical trajectories shown by dashed-dotted 
and dashed curves. We fitted the data by a pair of simple 
Gaussian functions and we obtained the dispersion rela- 
tions in Fig. [2] (b). Two curves from crystals A and B are 
separately observed, which are identical by translational 
displacement of h' = —{h + 0.53/). Thus the dispersion 
in the a* direction with the boundary energy of about 3 
meV is successfully obtained in Cu2Fe2Ge40i3. In com- 
bination with the previous datai^, the whole dispersions 
in three independent directions are drawn in Fig. [3] (a). 
Filled symbols are the data by single crystal NIS and 
the gray shaded area around 24 meV is a narrow band 
excitation observed in powder NIS experiment. 



In the 1st scenario we analyze the data by the con- 
ventional SW model. Since all the spins are confined in 
almost a - c plane and the canting is smalU^, we start 
from coUinear structure. We assume isotropic interac- 
tions for Fe - Fe, Cu - Fe, and Cu - Cu spins for the 
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TABLE I: Exchange parameters for Cu2Fe2Ge40i3 determined by neutron inelastic scattering (NIS) and neutron diffraction 
(ND). The values in the square bracket are the estimates by semi-static method. 





Jcu (meV) 


Jb (meV) 


Jc (meV) 


Jeff (meV) 


JCu-Fc/Jcu 


Jcu-Fo (meV) 


A (meV) 


NIS and ND (Refs. U and 17) 


2.4(2) X 10 


1.60(2) 


0.12(1) 


[0.13(4)] 


[0.10(5)] 


[2.4(2)] 




NIS (Conventional SW) 


2.(4) X 10^ 


1.59(2) 
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FIG. 3: (Color online) (a) The dispersion relation in 
Cu2Fe2Ge40i3. The lower energy excitation in o* direction 
(filled circles) is determined in the current study. Other data 
are from Ref. [itI . Solid lines and small open circles are the fit 
by conventional SW and effective SW model (see text), (b) 
Energy scales of Fe chains and Cu dimers. Dynamic suscep- 
tibility of dimers is indicated by a solid curve. 

simplicity. Hence the Hamiltonian is defined as 

1 / T oFc (-.Fo I T qFo ^Cu 

At = - 2^ yJaps^ ■ bp + j^s^j ■ J- S 

Q,/3,-ve8 Fe sublatticcs, 
5,e,Ce8 Cu sublattices 

+ J,c^f"-T^"). (1) 

Here a, /3, and 7 denote Fe sublattices and S, e, and ( de- 
note Cu sublattices. and T^" are the spin operators 
of S' = 5/2 and S = 1/2 on Fe and Cu ions, respec- 
tively. We considered Jb, Jc, Jcu-Fe, and Jcu in Fig. [1] 
(b) as the contributed exchange parameters in eq. ([T]). 
In addition the split anisotropy gaps A^; and Aj, are em- 
pirically introduced as (hjJx,y{<l))^ — (^(q))^ + 
^x,y =0.95 and 1.16 meV were obtained by fitting the 
resolution convoluted function to the constant q scan at 
g = (2 1 -0.5), where the xy anisotropy is the most en- 
hanced. The split is small enough and we will use the 
mean value A = 1.0(5) meV as fixed parameter here- 



after. 

The numerical fit to the dispersion is indicated by small 
open circles in Fig. [3] (a). The fit to the data is very good. 
While doing simulation it is found that dispersionless ex- 
citation at fiw 24 meV, where four modes are over- 
lapped, is sensitive to the change of parameter Jcu- This 
means that they are from Cu centered spins. Similarly 
four modes in the lower energy are from Fe centered spins. 
The estimated exchange parameters relating to Fe spins 
are summarized in the second line in Table HI They are 
reasonably consistent with previously estimated values in 
the first line. However, Jcu of 2.(6) x 10^ meV is very 
different from the previous value 24.(2) meV in Ref. 
This is because of the different assumption for the ground 
state; classical Neel state for the present SW and singlet 
dimer state for the previous analysis. To discuss the Jcu 
the isostructural compound Cu2Sc2Ge40i3^1^ where Fe 
ions in Fig. [T] (b) are replaced by nonmagnetic Sc ions 
would give good insight. In these compounds the angle 
Z Cu-O-Cu and the bond distance are very close to each 
other and so would the magnitude of superexchange con- 
stant be. The intradimer interaction in Cu2Sc2Ge40i3 
was reliably determined as Jsc = 24.(5) ^ 25.(4) meV 
by the isolated dimers modeli^. The value is close to the 
previous estimate of Jcu 24.(2) meV— rather than the 
present one 2.(6) x 10^ meV. Hence, the singlet dimer 
ground state is a better assumption and we adopt the 
Jcu ~ 24.(2) meV as the reasonable value hereafter. 

In the 2nd scenario we will consider the effective SW 
model based on the energy scale separation of the mag- 
netic excitationsii. The absence of dispersion in the 
Cu centered mode and the suppressed magnetic moment 
means that Cu spins are strongly bound to spin dimers 
singlet state. Onto the Cul spin in the dimer the field, 
/iCui, is applied through Jcu-Fe by the Fel spins in Fig.[T] 
(c). Then the polarized spins, S'cu = x(7r,w)ft.cu, are in- 
duced on the both Cul and Cu2 ions that transfer the 
exchange integral from one Fe spin to another. Here the 
dynamic susceptibility of the dimer x(7r, lo) is purely real 
and expressed by x(7r,a;) = l/(Jcu ~ f^) at T = K 
(solid curve in Fig. [3] (b)). In the energy region of Fe 
spin networks (gray shaded area) x(7r, w) is almost 1/Jcu 
and constant. This suggests that static treatment of MF- 
RPA works well when we consider the spin dynamics at 
hoj < 10 meV. Then the indirect coupling Jeff can be ap- 
proximately expressed as Joff(ix') — J^^__pj,x(7r, a;)/2 ~ 
Jq^_p^/2Jcu- Thus the effective spin Hamiltonian for 
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the low energy region is 

^ = 2 X! JapS^J" ■ S^". (2) 

a,/3G8 Fc sublatticcs 

I 



H ere A = 2Jb + JcS + 2Jc, B = 
\J Jl ^ J^g + 2 Jfc Joff cos 2'Kh, C = Jb, and 
D = 2 Jc cos 27rL Reasonable fit to the data is shown by 
the thick solid curves in Fig. [3] (a). They are perfectly 
overlapped on the open small circles calculated in the 
1st scenario. The obtained exchange parameters are 
summarized in the third line in Table HI 

Excellent consistency was found in the independently 
estimated exchange parameters of Fe-centered spins. 
Particularly interesting is the effective indirect coupling 
Jeff's. In the current study we obtained Jeff = 0.09(6) 
meV by the direct measurement of the spin dynamics. 
In the previous study^^ Jcu-Fo/Jcu was estimated stat- 
ically in MF-RPA level from the staggered magnetiza- 
tion curves of Cu dimers obtained by ND. Combining 
the parameter with Jcu = 24.(2) meV in a separate 
experimentii, J^s = 0.13(4) meV was estimated semi- 
statically. The reasonable consistency in Jeff's means 
that the static MF-RPA treatment works well and, thus, 
Cu dimers behave as the nonmagnetic media in the indi- 
rect Fe-Fe coupling. Due to the simplicity of the model, 
all the related exchange constants, Jcu, Jeff, and Jcu~Fc 
are experimentally obtained. 

There are a few other compounds reported as bicom- 
ponent systems with spin-singlet entity and additional 
spin. Among them i?2BaNi05 (i? = rare earth metals) 
have been studied the most in detaiU^. The compound is 
described by reasonably decoupled one dimensional Ni^+ 
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We considered Jb, Jc, and JcS as the exchange parame- 
ters as is illustrated in Fig. [T] (e) . The calculated disper- 
sion relation is 



(3) 

I 

{S = 1) chains and surrounding rare earth ions. The 
weak coupling induce magnetically ordered state of which 
the nature is similar to that in Cu2Fe2Ge40i3^*'. How- 
ever, the rare earth moment behaves as Ising type spin 
and the indirect interaction by way of Haldane chains has 
not been studied. Other examples are Cu2CdB2 06 and 
Cu3(P2060H)2^^. Spin model of quantum dimers and 
nearly free spins is proposed from the observed magneti- 
zation plateaux in the M - H curves. For lack of single 
crystal, however, the detailed spin dynamics is still un- 
known. As far as we know Cu2Fe2Ge40i3 is the only 
model compound for the indirect magnetic interaction 
by way of spin singlet media. 

In conclusion we succeeded to obtain the spin dis- 
persion in the a* direction in Cu2Fe2Ge40i3 by high- 
resolution neutron scattering experiment and by the elab- 
orate data analysis. The whole dispersion is well ex- 
plained by the effective SW model. The effective Fe - Fe 
interaction Jeff is directly estimated and the value is con- 
sistent with that obtained semi-statically. It is found that 
Cu2Fe2Ge40i3 is a prototype compound for randomness- 
free RKKY-like interaction in insulating metal oxides. 
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